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ABSTRACT We have recently reported the isolation of a novel virus, provisionally designated C/swine/Oklahoma/1334/201 1 (C/ 
OK), with 50% overall homology to human influenza C viruses (ICV), from a pig in Oklahoma. Deep RNA sequencing of C/OK 
virus found a matrix 1 (M 1 ) protein expression strategy that differed from that of ICV. The novelty of C/OK virus prompted us 
to investigate whether C/OK virus could exist in a nonswine species. Significantly, we found that C/OK virus was widespread in 
U.S. bovine herds, as demonstrated by reverse transcription (RT)-PCR and serological assays. Genome sequencing of three bo- 
vine viruses isolated from two herds in different states further confirmed these findings. To determine whether swine/bovine 
C/OK viruses can undergo reassortment with human ICV, and to clarify the taxonomic status of C/OK, in vitro reassortment 
and serological typing by agar gel immunodiffusion (AGID) were conducted. In vitro reassortment using two human ICV and 
two swine and bovine C/OK viruses demonstrated that human ICV and C/OK viruses were unable to reassort and produce viable 
progeny. AntigenicaUy, no cross-recognition of detergent split virions was observed in AGID between human and nonhuman 
viruses by using polyclonal antibodies that were reactive to cognate antigens. Taken together, these results demonstrate that 
C/OK virus is genetically and antigenicaUy distinct from ICV. The classification of the new virus in a separate genus of the Ortho- 
myxoviridae family is proposed. The finding of C/OK virus in swine and bovine indicates that this new virus may spread and es- 
tablish infection in other mammals, including humans. 

IMPORTANCE Influenza C viruses (ICV) are common human pathogens, infecting most people during childhood and adoles- 
cence, and typically cause mUd respiratory symptoms. While ICV have been isolated from both pigs and dogs, humans are 
thought to be the natural viral reservoir. Previously, we characterized an ICV-Uke virus isolated from pigs exhibiting symptoms 
of influenza virus-like illness. Here, we show molecular and serological data demonstrating widespread circulation of simUar 
viruses in bovines. Deep RNA sequencing, phylogenetic analysis, and in vitro reassortment experiments demonstrate that ani- 
mal ICV-like viruses are geneticaUy distinct from human ICV. AntigenicaUy, we show that ICV-like viruses are not recognized by 
ICV antibodies. En masse, these results suggest that bovine influenza virus warrants classification as a new genus of influenza 
virus. The finding of this novel virus that can infect multiple mammalian species warrants further research into its role in human 
health. 
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The Orthomyxoviridae family contains three genera of influenza 
viruses, influenza A, B, and C viruses, that are classified ac- 
cording to antigenic differences between their nucleoprotein (NP) 
and matrix 1 (Ml) proteins (1-3). Such a classification is sup- 
ported by low intergenic (20 to 30%) and high intragenic (>85%) 
homology of their Ml and NP proteins. Additionally, viral reas- 
sortment has been used to classif)' influenza viruses, as by defini- 
tion, intragenic reassortment yields viable progeny (4). 

Influenza C virus (ICV) is a ubiquitous pathogen of humans. 
Compared to influenza A and B viruses, it causes a mQd respira- 
tory illness or remains an asymptomatic infection. Most humans 
are seropositive to ICV by adolescence, and ICV is most com- 



monly isolated from children less than 6 years of age (5, 6). Clin- 
ically, fever and respiratory symptoms are common and typically 
mild. Aside from humans, ICV has been isolated from pigs in 
China in 1981 and in the United States in 2011 (7, 8). Genetic 
analysis of the Chinese swine ICV found high similarities to con- 
temporary human ICV, suggestive of interspecies transmission 
(9). In contrast, the ICV-like C/swine/Oklahoma/1334/2011 (C/ 
OK) virus was only moderately related to previously characterized 
ICV, with approximately 50% overall homology. These results, in 
addition to a lack of cross-reactivity between antibodies against 
C/OK virus and human ICV in hemagglutination inhibition (HI) 
assays, suggested that C/OK virus conservatively represents a new 
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but distant subtype of ICV. Previously, only influenza A viruses 
were thought to possess multiple subtypes (10). Cocirculation of 
multiple influenza virus subtypes allows for antigenic shift, the 
process of rapid evolution mediated by reassortment of the 
genomic RNA segments of two parent viruses. 

Serological studies found that 19% and 9.9% of pigs in the 
United Kingdom and Japan, respectively, were seropositive to hu- 
man ICV (11, 12). Similarly, 9.5% of pigs had measurable titers in 
the HI assay to C/OK virus ( 8 ) . Together, these results suggest that 
pigs are frequently infected by ICV. Despite its ability to infect 
ferrets, a surrogate for human influenza virus studies, only 1.3% of 
human sera had a measurable titer to C/OK virus. Additionally, 
the relatively low geometric mean titer (GMT) in swine sera in the 
HI assay to C/OK virus (20.7) led us to speculate that neither swine 
nor humans are the natural reservoir of C/OK virus. 

Using a reverse transcription (RT)-PCR primer and probe set 
targeting the FBI gene of swine C/OK virus, we identified C/OK 
virus in samples from bovine respiratory disease diagnostic sub- 
missions. Detailed genetic, biological, and antigenic characteriza- 
tion suggests that C/OK viruses in cattle and swine are distinct 
from human ICV and likely represent a new genus in the Ortho- 
myxoviridae family, with C/OK virus as the type species of influ- 
enza D virus. Additionally, widespread and high antibody titers to 
C/OK virus and relative ease of virus isolation in cattle indicate 
that bovines may represent the reservoir to this novel influenza 
virus. The finding of C/OK virus in both swine and bovine indi- 
cates that this new virus may spread and establish infection in 
other mammals, including humans. 

RESULTS 

Analysis of viral transcriptome by next-generation sequencing 
(NGS). To extend the understanding of swine C/OK virus genome 
expression, we used lUumina-based deep RNA sequencing to 
study viral transcription in swine testicle (ST) cells, focusing on 
examining and comparing the transcription and splicing of swine 
C/OK virus with those reported previously for human ICV. Total 
cellular RNA was isolated from ST cells infected with swine C/OK 
at a multiphcity of infection (MOI) of 0.1 at 0, 18, and 36 h after 
infection. Corresponding cDNA libraries were prepared and 
loaded onto an lUumina Genome Analyzer IIx (GA IIx) sequencer 
for sequencing. The DNA sequence data sets were then assembled 
against the C/OK genome, which has seven single-stranded, 
negative-sense RNA segments. For each sample, we obtained on 
average over 9 million 100-nucleotide (nt) reads mapping to ei- 
ther C/OK virus or pig genomes. At 18 h postinfection, approxi- 
mately 5% of the total reads mapped to the viral genomic seg- 
ments, which was similar to that obtained at 36 h postinfection 
(data not shown). Read coverage for each of seven segments had 
significant overlap at both time points, demonstrating the active 
replication of the C/OK genome with a similar rate at both time 
points (Fig. 1). The pattern ofread coverage clearly discriminated 
seven segments into three distinctive groups (Fig. 1). Group 1 
included three viral polymerase genes (P3, PBl, and PB2 genes) 
that all had a large proportion of reads mapping to the 5' and 3' 
ends of segments, with a significant reduction in the number of 
reads mapping to the internal regions. The overrepresentation of 
sequencing reads in the 5' and 3' ends of the segments suggested 
the presence of internal deletions in copies of these polymerase 
genes. The subgenomic segments with the internal deletions ap- 
peared to range in size from 300 to 700 nucleotides, which have 



been previously found to be associated with defective interfering 
RNA in replicating influenza viruses, including human ICV (13- 
16). Group 2, represented by NP and hemagglutinin-esterase fu- 
sion (HEF) segments, had sequencing reads relatively evenly dis- 
tributed across the complete segment compared to group 1. 
Interestingly, a reduction of read coverage was noted in the 5'- 
terminal half of the NP segment. The coverage reduction may 
indicate the presence of heterogeneous species of the NP genes 
with some sequences deleted within this region, which is similar to 
those observed for polymerase genes (PBl, PB2, and P3 genes). 
Group 3 (M and NS segments) displayed a large reduction in the 
number of reads mapping in certain regions of the segments, in- 
dicative of the presence of internal introns. Subsequent analysis of 
segment-derived sequencing reads identified potential splice 
junctions and introns present in both segments (Fig. 1; see Ta- 
ble SI and S2 in the supplemental material). 

Follow-up RT-PCR and Sanger sequencing experiments con- 
firmed the splice junctions for both NS and M segments and their 
introns determined in NGS analysis (see Fig. S2 to S4 and Ta- 
bles SI and S2 in the supplemental material). Further analysis 
showed that both splicing events utilized a consensus dinucleotide 
motif GT/AG: GT present at the extreme 5' (donor) with AG 
detected at the extreme 3' (acceptor) ends of the intron. The splic- 
ing pattern for the C/OK NS segment was similar to that of human 
influenza C virus (Fig. 2A) (17). In contrast, M segment splicing 
encoding mRNA for the Ml protein was different from that of the 
human ICV Ml mRNA (Fig. 2B). Rather than introduction of a 
termination codon to the Ml open reading frame, splicing of the 
C/OK M segment added a second exon encoding a 4-residue pep- 
tide to the primary Ml transcript (Ml mRNA, nucleotides 29 to 
756 joined to 1050 to 1062). Variability in the length of the second 
exon of Ml observed between the two viruses adds to the novelty 
of the C/OK virus. In addition to NS2 and Ml expressed through 
splicing, NS and M segments of both human ICV and C/OK vi- 
ruses also encode NSl and M2 proteins, respectively, through the 
corresponding unspliced mRNA transcripts. These viruses likely 
utilize a similar mechanism to produce these proteins, although 
further experimental validation is required. 

Identification of swine C/OK virus in bovines. The novelty of 
C/OK virus demonstrated above and biological distinctness from 
human ICV reported in our previous work (8) prompted us to 
investigate whether C/OK virus could exist in a nonswine species. 
Nasal swabs were submitted to Newport Laboratories for diagnos- 
tic testing from cattle with respiratory disease. A total of 45 sam- 
ples were tested by relative real-time reverse transcription-PCR 
(rRT-PCR) for C/OK virus. Samples were collected in the late 
winter and spring of 2013 and originated from 6 different states. 
Eight samples (18%) were positive, all from Minnesota and Okla- 
homa. Virus isolation was attempted on PCR-positive samples 
using multiple cell lines, and viral replication was determined by 
rRT-PCR and hemagglutination (HA) assays. Virus was initially 
isolated from five different samples on human rectal tumor 
(HRT-18G) cells. The virus was noncytopathic but replicated to a 
high titer as determined by rRT-PCR (threshold cycle [Cj] values 
of -12 to 14) and hemagglutination assays (HA = 1,280 to 5,120 
hemagglutinating units [HAU]/ml). Following the initial isola- 
tion, viruses were passaged on ST cells and displayed typical influ- 
enza virus cytopathic effects. Four of the isolated viruses all orig- 
inated from the same herd in Minnesota. Two were randomly 
chosen for further study, C/bovine/Minnesota/628/2013 (C/628) 
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FIG 1 Transci'iptome maps of C/OK virus segments. Numbers of normalized sequence reads (100 nucleotides each read) were mapped to each nucleotide 
position across the respective segments of C/OK virus genome. The x axis represents the length of each segment, while the / axis represents the number of 
sequencing reads mapping to each genomic segment. Arrow signs indicate intron regions present in NS and M segments that were identified by next- generation 
sequencing (NGS) analysis. 



and C/bovine/Minnesota/729/2013 (C/729), along with the re- 
maining viral isolate from a bovine in Oklahoma, C/bovine/Okla- 
homa/660/2013 (C/660). 

Viral genome sequencing. The genome sequences of C/628, 
C/729, and C/660 were determined using an Ion Torrent se- 
quencer. Primers designed to amplify the whole genome of lAV 
were modified with the 5 ' and 3 ' noncoding regions of C/OK virus 
(8, 18). Using viral genomic RNA as a template, the seven- 
segmented genome was amplified from the three viruses using a 
high-fidelity polymerase, and the resulting amplicons were se- 



quenced. Sequences were assembled using C/OK virus as a tem- 
plate. As amplicons were produced using primers homologous to 
the 5' and 3' noncoding regions of C/OK virus, the sequence en- 
compassed the segment between these regions. The full-length 
segments were sequenced with the exception of 5' and 3' noncod- 
ing regions. All segments had greater than 96% identity to C/OK 
virus. Segments with the highest homology to C/OK virus were 
PBl (98.9 to 99.1% identity) and NS (98.8 to 99.2% identity), 
while the most divergent segments were HEF (96.7 to 99.0% iden- 
tity) and P42 (96.9 to 99.2% identity). 
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FIG 2 Splicing strategies of C/OK virus for the NS segment (A) and M segment (B). Panel A scliematically illustrates a splicing strategy of the C/OK virus NS 
segment to produce NS2 protein in comparison to its counterpart in human influenza C viruses (ICV), while panel B describes a splicing strategy of the M 
segment to produce Ml protein in comparison to human ICV's Ml protein synthesis. 



RNA sequencing was performed to determine the conserved 5' 
and 3' noncoding region sequences for C/628, C/729, and C/660. 
The terminal conserved noncoding regions for each segment were 
identical to those of C/OK (data not shown). 

Phylogenetic analysis. To further understand viral evolution, 
we performed phylogenetic analysis of bovine C/OK virus and 
representative influenza A, B, and C viruses as well as other mem- 
bers of the Orthomyxoviridae family, Thogoto virus and infectious 



salmon anemia virus (Fig. 3). AH seven segments were analyzed 
individually using maximum likelihood analysis. All segments for 
C/628, C/660, and C/729 were closely related to C/OK. C/OK, 
C/628, C/660, and C/729 clustered most closely to human ICV, 
suggesting that these viruses diverged from human ICV after the 
divergence from ancestral lAV and IBV. Despite reports of fre- 
quent reassortment between ICV (19), all segments of the bovine 
ICV were highly homologous to C/OK, and there was no evidence 
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FIG 3 Phylogenetic trees of the seven segments of bovine C/OK viruses. Maximum likelihood analysis in combination with 500 bootstrap replicates was used 
to derive trees based on the nucleotide sequences of the genomic segments. Bootstrap values are shown above and to the left of the major nodes. 
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FIG 4 Mean hemagglutination inliibition titers for eiglit bovine herds ( 1 1 to 27 animals per herd) to bovine C/660 and swine CI OK. Bovine sera from eight herds 
located in five different states were screened individually ( 1 1 to 27 animals per herd) in the HI assay using C/OK and C/660 as antigens. The x axis represents the 
sample size of each herd tested, while the y axis indicates the average HI titer and standard error from each herd. 



of reassortment with human ICV. These findings led us to hypoth- 
esize that the C/OK clade of viruses is unable to reassort with 
human ICV. 

In vitro reassortment. To test our hypothesis and determine 
whether C/OK and the bovine ICV viruses are able to reassort with 
human ICV, an in vitro reassortment experiment was performed. 
Elucidation of reassortment between C/OK and human ICV en- 
abled us determine whether this novel virus utilizes reassortment 
as a potential mechanism for viral evolution and production of 
new antigenic variants that may pose a threat to human and ani- 
mal health. Another purpose of this experiment was to gain more 
evidence toward understanding the taxonomic status of C/OK 
virus, because, by definition, all members of an influenza virus 
genus can reassort with each other and give viable progeny (4). 

Two representative human ICV (C/Taylor/1947 [C/Tay] and 
C/Johannesburg/66 [C/JHB]), the swine C/OK, and the bovine 



C/660 were chosen. Selection of bovine C/660 over other bovine 
C/OK isolates is due to considerations that this virus encodes a 
slightly divergent HEF protein. ST cells were coinfected with a 
multiplicity of infection of 0.1 for each virus with all pairwise 
combinations of virus, with all four viruses, or, as a control, cells 
were also inoculated singly with each virus (Table 1 ) . Ten plaques 
from each well containing multiple viruses were analyzed by full 
genome sequencing to allow determination of segment parentage. 
The donor virus for each segment was determined for each 
plaque-purified virus using criteria of at least 99.5% identity. Due 
to the high similarity between the two human and also the two 
nonhuman viruses, the parent donor virus for every segment 
could not be determined. 

Reassortant viruses were identified following coinfection with 
the two human ICV, C/JHB and C/Tay, despite our inability to 
assign a parent virus for five of the 7 segments due to high se- 



TABLE 1 Genotypes identified in plaque-purified viruses isolated following infection of cells with influenza C viruses 



Parentage (no. of plaques)' 



Viruses" 


PB2 


FBI 


P3 


HEF 


NP 


P42 


NS 


C/IHB, C/Tay, 


C/OK (8); 


C/660 (9); 


C/OK (6); 


C/OK (7); C/660 (1); 


C/OK (4); C/660 (4); 


C/OK (1); C/660 (4); 


C/OK (7); C/660 (3) 


C/OK, C/660 


C/660 (2) 


NI(1) 


C/660 (4) 


NI(2) 


NI(2) 


NI(5) 




C/IHB, C/Tay 


NI (10) 


NI (10) 


C/JHB (7); 


C/Tay (4); NI (6) 


NI (10) 


C/Tay (10) 


NI (10) 








C/Tay (3) 










C/JHB, C/OK 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/JHB, C/660 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 ( 10) 


C/660 (10) 


C/Tay, C/OK 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/OK (10) 


C/Tay, C/660 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/660 (10) 


C/OK, C/660 


C/OK (10) 


C/OK (3); 


C/OK (8); 


C/OK (10) 


C/OK (6); 


C/OK (4); C/660 (5); 


C/OK (6); C/660 (1); 






C/660 (7) 


C/660 (2) 




C/660 (4) 


NI(1) 


NI (3) 



" Viruses used for coinfection. 

^ Parentage of viral genome segments present in virus plaques from coinfected cells. Ten plaque-purified viruses were analyzed from each coinfection experiment. The number of 
plaques from each donor is indicated in parentheses. Some viral segment donors could not be identified (NI). 
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TABLE 2 Results of AGID assay 

Presence or absence of visible white precipitation line" 
Antigen (virus or A/NWS/34(H1)- B/Hong 

mock control) A/Equine/Prague/ 1/56(N7) Kong/8/73(Matrk) C/Taylor/1233/1947 C/swine/OK/ 1334/20 11 

A/WSN/1933 + _ _ _ 

B/Brisbane/60/2008 - + - - 

C/Taylor/1233/1947 - - + - 

C/Johannesburg/1/1966 - - + - 

C/swine/OK/ 1334/20 11 - - - + 

C/bovine/OK/660/2013 - - - + 

MDCKmock* - _ _ _ 

HRT-18Gmock - _ _ _ 

PBS - _ _ _ 

" +, presence of a visible white precipitation line between antigen and antiserum wells; — , absence of a visible white precipitation line between antigen and antiserum weUs. 
The same protocol to produce viral antigens was employed to prepare supernatants collected from uninfected cells. 



quence identity. Only the P3 and P42 segments were divergent 
enough to conclusively determine lineage. All P42 segments were 
derived from C/Tay, while only 3/10 viruses had P3 segments de- 
rived from C/Tay. Reassortant viruses were also identified follow- 
ing coinfection of cells with the two nonhuman ICV, C/OK and 
C/660, representing at least 7 distinct genotypes (see Fig. SI A in 
the supplemental material). In contrast, coinfection of cells with 
pairwise combinations of human and nonhuman ICV failed to 
yield reassortant viruses. In all four instances, all segments from all 
viruses were derived from the nonhuman ICV (C/OK or C/660). 
Consistent with these results, coinfection of cells with all four 
viruses gave reassortant viruses with segments derived from C/OK 
and C/660 (see Fig. SIB). 

A second reassortment experiment was performed with C/JHB 
and C/OK using increased MOIs and an increased ratio of C/JHB 
to C/OK. Cells were coinfected with both C/JHB and C/OK with 
an MOI of 1.0:0.1, 1.0:1.0, or 10.0:1.0; as in the first reassortment 
experiment, all segments from all plaques isolated from cells coin- 
fected with C/JHB and C/OK were derived from C/OK. Eighteen, 
16, and 18 plaques, respectively, were analyzed by genome se- 
quencing. Similar to the previous reassortment experiment, all 
segments were derived from C/OK. 

In vitro growth kinetics were determined on ST cells for C/OK 
and C/JHB for the MOIs used in the second reassortment experi- 
ment (see Fig. S5 in the supplemental material). Cell inoculation 
was performed as in the reassortment experiment. Both C/JHB 
and C/OK had similar growth curves, with C/JHB consistently 
replicating at higher titers, approximately 1.0 log^ 50% tissue 
culture infectious dose (TCID5o)/ml, than C/OK at 24, 60, and 
72 h postinfection with the MOIs used in the second reassortment 
experiment. Despite these results, all segments from the second 
reassortment experiment were derived from C/OK, and no reas- 
sortant viruses were identified. 

The finding of reassortant viruses between the human viruses 
C/JHB and C/Tay as well as between the animal viruses C/OK and 
C/660 and the lack of reassortment between human and animal 
ICV suggests that human and animal ICV are unable to reassort 
and generate viable progeny. The failure to identify viable reas- 
sortant viruses provides the first functional evidence that this new 
group of viruses is genetically distinct from human ICV and rep- 
resents a new genus of influenza virus in the Orthomyxoviridae 
family. 

Antibody cross-reactivity. The inability of swine C/OK and 
bovine C/660 to reassort in vitro with human ICV plus the ob- 



served novelty of C/OK virus prompted us to employ serological 
typing by agar gel immunodiffusion (AGID) to probe their anti- 
genic relatedness. Antigenic differences between NP and Ml pro- 
teins detected in AGID or other similar assays allow division of 
influenza viruses into three genera (3, 20). AGID has been widely 
and routinely used to type influenza viruses, especially prior to the 
advent of RT-PCR and sequencing technologies. For the AGID 
experiment, similar amounts of split viral antigens (12 logj HAU) 
were used with normal uninfected cell-derived mock antigens as 
negative controls. Antisera generated against representative lAV, 
IBV, and ICV recognized intragenus, not intergenus, virus anti- 
gens in the AGID (Table 2). No cross-reactivity between genera 
and no recognition of mock antigens were observed, indicating 
the specificity of the AGID. Significantly, polyclonal antibody 
against C/Tay virus did not recognize detergent-treated swine 
C/OK and bovine C/660 virus antigens but was reactive to C/JHB 
virus antigen, in addition to its own antigen. Conversely, both 
swine C/OK and bovine influenza C virus antigens, not human 
influenza C virus antigens (C/JHB and C/Taylor) , were recognized 
by C/OK virus-specific polyclonal antibody. Considering that NP 
and Ml are major structural proteins of influenza virus, we con- 
cluded that the major differences of NP and M 1 are responsible for 
the failure in cross-recognition of viral antigens between these two 
groups of viruses in the AGID. This result is similar to what we 
observed for reference influenza A, B, and C viruses in the AGID. 
As such, based on this finding together with results from in vitro 
reassortment and NGS analysis, we propose that C/OK virus may 
be considered a representative strain of a new genus in the Ortho- 
myxoviridae family. 

Seroprevalence of C/OK virus in bovines. Given the relatively 
high percentage of bovine respiratory disease samples positive for 
C/OK-like virus by rRT-PCR, HI assays were performed on bo- 
vine sera to determine if antibodies to this virus were similarly 
prevalent. Our previous study of seroprevalence of C/OK virus in 
swine found 9.5% of pigs were positive, with a GMT of 20.7 (8). 
Bovine sera from eight herds located in five different states were 
screened individually (11 to 27 animals per herd) in the HI assay 
using C/OK and C/660 as antigens (Fig. 4). With the exception of 
one herd, all herds had geometric mean titers greater than 40 to 
both viruses, suggesting that C/OK virus infection is widespread in 
bovines and that antibodies generated against C/OK and C/660 
are cross-reactive in the HI assay. The widespread presence of 
C/OK virus in bovines suggests that they likely represent a reser- 
voir of this new group of viruses. 
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DISCUSSION 

Until recently, ICV was thought to consist of a single subtype and 
possess the lowest rate of evolution among the various genera of 
influenza viruses (10). ICV, unlike lAV, is thought to exist in a 
human reservoir and has previously only rarely been shown to 
infect alternative species, such as pigs and dogs (7, 9, 21). Despite 
the presence of multiple cocirculating lineages, ICV are relatively 
genetically homogenous and thought to be near evolutionary 
equilibrium in humans (22). A recently isolated influenza C-like 
virus raised many questions about whether our understanding of 
ICV was correct. This virus, C/OK, was isolated from a pig in 
Oklahoma displaying influenza virus-like illness (8). Genetic anal- 
ysis found only moderate (-53%) similarityto human influenzae 
viruses. Additionally, serological analysis found that antibodies 
against C/OK failed to recognize lAV, IBV, or human ICV. This 
virus also displayed phenotypes of broad cell tropism and had 
permissive growth at 37°C, which further distinguish this novel 
virus from human ICV. C/OK was provisionally identified as a 
subtype of ICV primarily due to its overall genetic organization (7 
RNA segments) and both coding and noncoding region similari- 
ties to ICV. Despite these, the numerous incongruities observed 
with C/OK warranted further investigation. 

Previous serological results using both human and swine sera 
found that approximately 1 and 10%, respectively, were positive 
for antibodies to C/OK, with low titers (HI range of 10 to 40). 
Infrequent isolation of C/OK virus from pigs plus a lower percent- 
age of RT-PCR samples positive to C/OK virus (<0.1%) during 
the routine testing of nasal swabs from pigs challenged a theory 
that that swine is the primary host of C/OK virus. Based on these 
results, we suspected an alternate reservoir for this virus. Screen- 
ing of a bovine serum panel representing animals from multiple 
states and production conditions (cow-calf and dairy) identified 
widespread, high antibody titers to C/OK viruses, suggesting that 
bovines may represent the reservoir to this virus. Serologically, the 
magnitude and percent positive for HI titers were similar to values 
reported for wild waterfowl, the established reservoir of lAV (23). 
Previous studies have shown that bovines are susceptible to lAV 
infection, and multiple subtypes of lAV have been isolated from 
cattle (24). Correlations have also been shown between lAV anti- 
bodies and respiratory disease and reduced milk yield (25, 26). 
Antibodies to lAV, IBV, and ICV have all been reported in cattle 
(27); however, the percentage of infected animals and titers were 
considerably lower than reported here for C/OK virus. 

ICV differs from lAV and IBV in that it codes for a hemagglu- 
tinin esterase fusion protein that catalyzes receptor binding, cleav- 
age, and membrane fusion. As opposed to lAV and IBV, which use 
a-2,3- and a-2,6-linked sialic acids on cell membrane proteins as 
receptors, ICV binds to 9-0-acetyl-N-acetylneuraminic acid 
(Neu5,9Ac2) as a receptor (28, 29). Another bovine pathogen, 
bovine respiratory coronavirus (BCV), encodes an evolutionarily 
related hemagglutinin esterase that also utilizes Neu5,9Ac2 as a 
receptor (30, 31). BCV can infect both the upper and lower respi- 
ratory tract and cause respiratory disease, implying that the recep- 
tors required for ICV are present in the bovine respiratory tract 
(31). Pathogenesis studies are required to evaluate whether C/OK 
virus utilizes Neu5,9Ac2 as a functional receptor and can replicate 
throughout the respiratory tract and establish whether C/OK virus 
causes disease in cattle. 

The genome sequences of three bovine viruses were very sim- 



ilar to those of the recently described swine C/OK virus (8). Phy- 
logenetic analysis showed that this clade of viruses was most 
closely related to human ICV; however, the distance between an- 
imal C/OK viruses and human ICV was similar to the differences 
between lAV and IBV for many segments. Additionally, there was 
no evidence of reassortment between human ICV and nonhuman 
C/OK virus. While our sample set was limited to three bovine 
viruses, we suspected that C/OK virus might be unable to reassort 
with human ICV. The abilities to reassort and yield viable progeny 
are necessary criteria for inclusion of a virus within a genus of 
influenza virus (4). An in vitro experiment was designed to deter- 
mine if C/OK virus could generate viable reassortants with human 
ICV. Not unexpectedly, no reassortment between swine/bovine 
C/OK viruses and human ICV was observed despite the method- 
ologies' success in detecting reassortment separately between hu- 
man and nonhuman ICV. These results strongly suggest that 
C/OK viruses are unable to reassort with human ICV and by def- 
inition that C/OK viruses represent a new genus of influenza virus. 
Previous work found that the 5' and 3' noncoding regions of 
swine C/OK virus were similar to those of human ICV but did 
contain one difference (position 5 from the 3' terminus) and a 
polymorphism at position 1 from the 3' terminus (8). The con- 
served 5' and 3' noncoding regions for the three bovine viruses 
presented in this work were identical to those of C/OK. The non- 
coding regions are highly conserved within each influenza virus 
genus and form panhandle structures due to inverted comple- 
mentarity, which is critical for genome replication and packaging 
(32-35). The inability of C/OK viruses to reassort with human 
ICV maybe due to these minor differences in noncoding regions. 

Influenza viruses have long been classified into three genera on 
the basis of antigenic differences among the conserved NP and M 
proteins (4, 20, 36). Despite molecular tools such as RT-PCR and 
sequencing, the serological assay is stUl an important and irre- 
placeable method to type influenza viruses, especially for newly 
identified influenza viruses, such as C/OK virus in this case. Ho- 
mologies of NP and Ml between lAV and IBV are approximately 
30 to 40%, which are slightly higher than those observed between 
lAV and ICV and between IBV and ICV ( 1 0 to 20%) . The approx- 
imate 40% similarity between C/OK and human ICV Ml and NP 
proteins further confounds the taxonomic position of C/OK virus. 
Results of our AGID experiments suggested that the new group of 
viruses (swine C/OK and bovine C/660) is antigenically distinct 
from human ICV, which is in good agreement with results ob- 
tained from the in vitro reassortment experiment. 

Influenza viruses exploit various splicing strategies to express 
their viral proteins to maximize the genome coding capacity. To 
date, splicing has been identified only in the two smallest seg- 
ments, M and NS. It is generally believed that all influenza viruses 
possess a similar mechanism to generate spliced mRNA from the 
NS segment to code for the NS2 protein (17, 37). Different mech- 
anisms for M segment splicing toward the production of Ml or 
M2 protein are reported in lAV and ICV, while the M segment of 
IBV is not spliced. lAV utilizes an alternative splicing mechanism 
to synthesize its M2 protein, while translation of Ml protein is 
made from an unspliced colinear mRNA (38). In contrast, the Ml 
protein of ICV is produced from a spliced mRNA, while its M2 
protein is derived from a protein (p42) generated from an un- 
spliced colinear mRNA followed by proteolytic cleavage by a sig- 
nal peptidase at an internal signal sequence (2, 39, 40). Using deep 
RNA sequencing coupled with conventional RT-PCR and Sanger 
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sequencing for confirmation, alternatively spliced mRNA for both 
NS and M segments in C/OK-infected cells was identified. C/OK 
virus is similar to other influenza viruses in NS segment splicing to 
generate the NS2 protein. In contrast to the ICV M segment that 
produces the Ml protein through splicing that solely introduces a 
termination codon (consisting of nucleotides 752, 753, and 982) 
into a spliced mRNA, splicing of the C/OK M segment adds an 
additional 4-amino-acid peptide into the preceding exon. Vari- 
ability in the length of the second exon of Ml adds to the novelty 
of this newly identified virus in swine and cattle. 

It is unknown if C/OK virus has an impact on human health. 
The ability to infect and transmit in ferrets, a model for human 
influenza virus pathogenesis studies, suggests this pathogen has 
the potential to cause disease in humans. Limited serology on a 
subset of human serum samples showing a 1.3% positive rate also 
supports this hypothesis. The prevalence of C/OK virus in cattle 
and presumptive spillover to swine, both of which live in close 
proximity to humans, further highlights its potential threats to 
human health, which merit further studies. 

In summary, several lines of evidence reported here present a 
compelling argument that C/OK virus is sufficiently distinct from 
human ICV to warrant classification as a new genus of influenza 
virus. Phylogenetic analysis found that the divergence between 
swine/bovine C/OK viruses and human ICV is similar to that be- 
tween lAV and IBV. Serologically, C/OK virus is antigenically un- 
related to human ICV. C/OK virus also encodes a novel mecha- 
nism for generating the Ml protein and, importantly, is unable to 
reassert with human ICV and generate viable progeny. In ecology, 
this new group of viruses differs from human ICV in that it resides 
in bovine. Furthermore, a serological survey suggested that C/OK 
virus frequently replicates and establishes its infection in bovines, 
with both percent positive and geometric mean titers considerably 
higher than those observed for humans and pigs, suggestive that 
bovines represent the reservoir of C/OK-like viruses. In conclu- 
sion, we propose that C/OK virus warrants creation of a new in- 
fluenza virus genus, influenza D virus, with D/OK (formerly 
C/OK) as the type species. 

MATERIALS AND METHODS 

Virus isolation. Nasal swabs or lung tissues were collected from cattle 
exhibiting respiratory disease and submitted to Newport Laboratories for 
diagnostic investigations during the late winter and spring of 2013. Sam- 
ples originated from Oklahoma, Minnesota, Iowa, Indiana, Texas, Ne- 
braska, and Tennessee and originated from animals of different ages. Vi- 
rus isolation was performed on both ST and HRT- 1 8G cells. Virus growth 
was determined by cytopathic effects, in addition to rRT-PCR and hem- 
agglutination assays, using turkey red blood cells. rRT-PCR was per- 
formed as previously described (8). 

Genome sequencing and analysis. Viral RNA was isolated using a 
MagMAX viral RNA isolation kit. Full-genome amplification was per- 
formed as previously described except that the primers were modified to 
match the noncoding regions of C/OK: 3' end, 5' ACGCGTGATCGCAT 
AAGCAG 3,' and 5' end, 5' ACGCGTGATCAGCAGTAGCAAG 3' (18). 
Amplicons were used for library preparation using the NEBNext Fast 
DNA library prep set for Ion Torrent. Libraries were sequenced using an 
Ion Torrent personal genome machine with the manufacturer's reagents 
and protocols. Contigs were assembled using C/OK (accession numbers 
NC_922305 to NC_922311) and C/Ann Arbor/1/50 (accession numbers 
NC_006306 to NC_006312) as the templates and utilizing the SeqMan 
NGen module from DNAStar. Phylogenetic analysis was performed using 
Mega 5. 1 software using maximum likelihood analysis with 500 bootstrap 
replicates to verify tree topology. The terminal segment sequences of 



C/660, C/628, and C/729 were determined using an RNA sequencing kit 
(Ion Total RNA-Seq kit v2) by following the manufacturer's instructions 
with viral RNA as described above. 

Next-generation sequencing of the viral transcriptome. Swine testi- 
cle cells were infected with swine C/OK at an MOI of 0. 1 . Total cellular 
RNA was isolated from mock-infected cells (0 h) and infected cells at 18 h 
and 36 h postinfection, respectively, with TRIZOL followed by RNA pu- 
rification with ethanol. RNA quality was tested using a Bioanalyzer 2100 
(Agilent, Palo Alto, CA) and then processed for cDNA library construc- 
tion by using a cDNA library prep kit (catalog number FC-122-1001; 
lUumina Inc.) according to the manufacturer's instruction. cDNA library 
construction and all other procedures were conducted in Genomics Core 
Research Facility (GCRF) at the University of Nebraska- Lincoln (UNL). 
Briefly, mRNAs were purified from the total RNA using oligo(dT) mag- 
netic beads followed by fragmentation. The resultant mRNAs were reverse 
transcribed to cDNAs that were subjected to an end repair process fol- 
lowed by ligation to the adapters. After separation in agarose gel through 
electrophoresis, cDNA fragments with a size of about 200 bp were excised, 
extracted, and amplified by PGR using two primers that match the ends of 
adaptors. PCR-enriched samples were then sequenced by an lUumina 
Genome Analyzer IIx (GA IIx) sequencer in the GGRF at UNL. 

An infection experiment consisting of three time points (0, 1 8 h, and 
36 h) was duplicated, and derived RNA samples were sequenced in deep 
RNA sequencing analysis. Three samples representative of three time 
points from each experiment were combined and loaded into one lane for 
sequencing. Each of the 6 samples submitted had approximately 9 million 
100-nucleotide single-end reads, and the sequencing results passed qual- 
ity control. 

Bioinformatic analysis of the NGS data set. Reads from each dupli- 
cate were merged and mapped to C/OK reference genomes using the 
GSNAP program (41). The resulting SAM files were converted to BAM 
files, and read coverage at each position of the G/OK genome segments 
was calculated using SAMtools (42). Sequencing reads were normalized 
according to total number of reads to viral and pig genomes obtained at 
each time point. Integrative Genomics Viewer (IGV) was used for visual- 
ization and mapping quality assurance (43). Reads mapped to G/OK ge- 
nome segments with splicing were further analyzed with bash scripts. The 
accession numbers for the SRA database of NCBI are as follows: 
SRX471882 for the two repeats of 0 h, SRX471884 for 18 h, and 
SRX471846 for 36 h following infection. 

Cell culture and in vitro reassortment. ST and human rectal tumor 
(HRT-18G) cells were grown in Dulbecco's modified Eagle's medium 
(DMEM) containing 5% fetal bovine serum. Influenza G/Taylor/1233/ 
1947 (C/Tay) virus was provided by BEI Resources, while G/Iohannes- 
burg/1/66 (C/IHB) virus was a gift from P. Palese (Mt. Sinai Medical 
School, New York) . For routine cultivation, cell medium was replaced by 
DMEM and virus was inoculated with a multiplicity of infection (MOI) of 
0.00 1 . Nonhuman IGV were cultivated at 33''G or 37°G, while human IGV 
were cultivated at 33°G. For in vitro reassortment, all viruses were propa- 
gated at 33°C in ST cells and titrated to determine 50% tissue culture 
infectious dose per milliliter (TGID5o/ml). A single well of a 6-well plate 
with a monolayer of ST cells was coinfected with viruses as shown in 
Table 1 with an MOI of 0.1. Viruses were adsorbed for 1 h, and then the 
monolayer was washed with phosphate-buffered saline (PBS) and refed 
with DMEM. After 3 days, cell culture supernatants were serially diluted 
and plated onto a monolayer of ST cells with an agar overlay, and after an 
additional 1 to 3 days, individual plaques were picked and expanded on ST 
cells. Ten plaques were expanded for each reassortment experiment, and a 
single plaque was chosen from cultures of parental viruses. 

A second reassortment experiment was performed under conditions 
designed to optimize reassortment between G/IHB and G/OK. Cells were 
coinfected with both C/IHB and G/OK with either an MOI of 1.0:0.1, 
1.0:1.0, or 10.0:1.0, as described above. After 3 days, cell culture superna- 
tants were serially diluted and plated onto a monolayer of ST cells with an 
agar overlay, and after an additional 1 to 3 days, individual plaques were 
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picked and expanded on ST cells. Twenty plaques were expanded for each 
reassortment experiment. 

In vitro growth kinetics were determined on ST cells for C/JHB and 
C/OK with the MOIs used in the second reassortment experiment. Cell 
inoculation was performed as described above using the same viral stocks 
used in the reassortment experiments. Samples were collected at 24, 60, 
and 72 h postinfection and analyzed by titration on ST cells. 

Serology. Bovine sera were treated with receptor-destroying enzyme 
overnight at 37°C and then heat inactivated at 56°C for 30 min. Hemag- 
glutination inhibition (HI) assays were run with 0.5% packed turkey red 
blood cells as described in the WHO's Manual on Animal Influenza Diag- 
nosis and Surveillance (36). Bovine sera were a gift from the SD State 
University Veterinary Diagnostic Laboratory and represented samples 
from five states (South Dakota, Vermont, Pennsylvania, Idaho, and Cali- 
fornia) submitted for unrelated diagnostic testing. 

All influenza A, B, and C reference virus antigens and antisera were 
obtained from Biodefense and Emerging Infections Research Resources 
Repository (BEI Resources). AGID antigens for swine C/OK virus, bovine 
influenza C-like virus, human influenza virus C/JHB, and human influ- 
enza virus C/Taylor were prepared according to the protocol of WHO and 
the World Organization for Animal Health (OIE). Supernatants from 
uninfected cell cultures (MDCK and HRT-18G) were prepared as mock 
controls using the same procedures as those used for production of viral 
antigens. Rabbit polyclonal antibody against swine C/OK virus was pre- 
pared using purified C/OK virus antigens. AGID plates were prepared 
according to the standard procedure, and assays were performed in trip- 
licate according to WHO/OIE protocol (36). Similar concentrations of 
viral antigens for all tested strains (approximately 12 logj HA units per 
0.025 ml) were used in parallel with mock controls for antigen recognition 
by antisera in AGID with two different dilutions (neat and 1:2). 

Nucleotide sequence accession numbers. Genome sequences were 
entered into GenBank under accession numbers KF425652 to KF425658 
(C/628), KF425659 to KF425665 (C/660), and KF425666 to KF425672 
(C/729). 
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